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Measuring Top Quark Polarization in Top Pair plus Missing Energy Events 
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The polarization of a top quark can be sensitive to new physics beyond the standard model. 
Since the charged lepton from top quark decay is maximally correlated with the top quark spin, 
it is common to measure the polarization from the distribution in the angle between the charged 
lepton and the top quark directions. We propose a novel method based on the charged lepton 
energy fraction and illustrate the method with a detailed simulation of top quark pairs produced in 
supersymmetric top squark pair production. We show that the lepton energy ratio distribution that 
we define is very sensitive to the top quark polarization but insensitive to the precise measurement 
of the top quark energy. 



Introduction: Events with a top quark pair plus miss- 
ing energy (tt+ Et) are promising channels in which 
to investigate models of new physics (NP) beyond the 
standard model (SM). Missing energy originates typically 
from non-interacting or otherwise invisible dark matter 
(DM) candidates in the NP models, along with neutrinos 
from SM decays. In these events the polarization of the 
top quark is sensitive to the chirality structure of the top 
quark's interaction with a postulated parent new heavy 
resonance and the DM. The top quark polarization might 
provide a new way to gain insight into NP models. Mea- 
surements of the top quark polarization tend to rely on 
the predicted angular correlation of the momentum of a 
charged lepton (from the top quark decay) with the top 
quark spin [lj]. However, this measurement is difficult in 
tt+ Et events because it is generally not possible to re- 
construct the top quark kinematics, i.e., to disentangle 
the kinematic effects of the DM particles from neutrinos 
that accompany the charged leptons in the top quark 
leptonic decay. 

In this Letter we define and examine the energy frac- 
tion of the charged lepton from the top quark as a novel 
measure of top quark polarization, without the require- 
ment of top quark reconstruction and knowledge of the 
dark matter mass and spin. We emphasize a few advan- 
tages of our energy ratio variable: (i) it is sensitive to the 
top-quark polarization; (ii) it not sensitive to the mass 
splitting between a heavy resonance parent and the DM 
candidate, provided that this splitting is not too small; 
(iii) the difference between the left-handed top-quark (tz,) 
and the right-handed (tn) is not sensitive to the spin of 
a heavy parent resonance or to the collider energy. 

We illustrate our method with top squark (t) pair pro- 
duction in the minimal supersymmetric extension of the 
SM (MSSM), pp -> ttX -> ttxxX, where \ denotes a 
neutralino (the DM candidate) . Once a signal for a t 
or another such NP candidate has been established, an 
important next step would be to determine the charac- 
teristics of its interaction with the SM particles. We do 
a full model simulation of pp — > ttX — > tixxX at the 



Large Hadron Collider energy 8 TeV, including typical 
experimental selection cuts |2j. From these events, we 
compute the energy fraction of the charged lepton from 
the top quark decays, show explicitly the relationship of 
this energy fraction to the top quark spin, and demon- 
strate what one may conclude about the nature of the 
interaction of the top quark and top squark from such 
data. In the MSSM, the i-polarization probes the t-t-x 
interaction and in turn the top squark mixing Q . 
The method: In the leptonic decay of a top quark, 
t —> bW + —> M + v, the correlation of the momentum 
of the charged lepton £ + with the polarization §t of the 
top quark, viewed in the top quark rest frame, takes the 
form (1 + s t z)/2, where z = cos 9 is the cosine of the 
angle between the top quark spin axis and the lepton 
momentum. For a boosted top-quark with energy E tl 
the distribution in the energy fraction xi = 2Eg/E t of 
the charged lepton becomes 
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Here, A = Ty/ /mw is the ratio of the VF-boson width 
and VK-boson mass, B = m^/m 2 ps 0.216 is the ra- 
tio of the W^-boson mass and top quark mass, and the 
limits of integration are z m - m = max[(l — l/j 2 x)//3, —1], 
z max = min[(l - B/j 2 x)/f3, 1] with 7 = E t /m t and 
(3 = \J\ — l/"/ 2 . The function Arctan is defined as 
arctan(ir) for x > while 7r + arctan(x) for x < 0. 

Figure 0Ja) displays the normalized energy fraction 
of the charged lepton from Eq. [T] for three top quark 
energies, E t = (250 ,500, 1000) GeV, for both left- 
handed and right-handed top quark decay. The im- 
portant point .J, [f| is that right-handed top-quarks ta 
(dashed curves) produce more energetic leptons than left- 
handed top-quarks th (solid curves), with the difference 
becoming more pronounced with increasing E t . 

We exploit the different dependence of tz and tjj, on 
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FIG. 1: (a)Distributions of the energy fraction xe 

of a charged lepton from top quark decay for E± = 
250, 500, 1000 GeV; (b) The ratio Rasa function of the 
cut threshold x c for E t = 250, 500, 1000 GeV. The solid 
lines represent left-handed top quark decay while the dashed 
lines represent right-handed top quarks. 



xi shown in Fig [3a) to measure the top quark polar- 
ization Q. We introduce a ratio 1Z as a quantitative 
measure of the energy fractions of ti, and tR, 
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This ratio is a function of the the cut threshold x c of 
the energy fraction xe.. We plot the lZ(x c ) distribu- 
tion in Fig. QJb) f° r three top quark energies, E t = 
(250 ,500, 1000) GeV, for both t L (solid curves) and t R 
(dashed curves). While the energy distribution Fig.QJa) 
varies with the top-quark energy, the TZ(x c ) distribution 
in Fig. [TJb) shows much less dependence. 

An analytic expression can be derived for lZ(x c ) in the 
limit j3 — > 1. It takes the form 
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for x c € (25, 2), where A t = (— for (t L ,t R ), re- 
spectively. For small x c , these expressions show that 
lZ(x c ) grows linearly with x c for t^, whereas lZ(x c ) grows 
as x\ for tji (A t = 1). 

The analytic expressions Eqs. [3] and [4] also explain 
why the curves for E t = 500 GeV (0 = 0.94) and 
E t = 1000 GeV (P = 0.99) almost overlap. For an en- 
ergetic top-quark, an important consequence is that the 
difference between lZ(x c ) for and t R is not sensitive 
to E t , i. e., the mass splitting between the parent heavy 
resonance and the DM candidate, as long as the mass 
splitting is not too small. The and t R curves are in- 
sensitive to the origin of the top quark in the collision, 
whether from a heavy fermion decay or from a scalar 
decay. In other words, 7Z(x c ) quantifies the top quark 
polarization but not the top quark origin. Moreover, in 
order to extract NP signal events from SM backgrounds, 
one must normally impose a set of hard kinematic cuts 
on the leptons and jets in the final state. These hard 
cuts force the top quark to be very energetic and thus to 
satisfy the limit @ — > 1. Therefore, another virtue of the 
7i(x c ) variable is that the difference between and t R 
curves, do not vary with the hard cuts. 

The ratio lZ(x c ) appears to show great promise for dis- 
tinguishing and ifl.. However, even if it is insensitive 
to Et, it presupposes reconstruction of the kinematics 
of the top quark (i.e., knowledge of E t ). Moreover, un- 
til this point, we have not included the influence of the 
production dynamics of the top quark, including matrix 
elements and the convolution with parton distribution 
functions. To prove our method useful, we must show 
that there are good estimators that can replace E t . To 
this end, we turn to an explicit calculation of top squark 
(t) pair production, pp —> ttX — > ti\xX. 
Collider simulation: We perform a parton-level Monte 
Carlo simulation of top squark (t) pair production pp — > 

ttX -+ ttxxX to demonstrate that 1Z remains useful for 
distinguishing tL and tR even when Et cannot be mea- 
sured directly. We assume the colored scalar t decays 
entirely into tx through the effective coupling 



g cS tx(cos cS P L + sin6 cS P R )t, 



(5) 



where the angle Q c s depends on the mass matrix mixings 
of the top squark and the neutralino sectors, and Pl/r 
is the usual left/right- handed projector. Our benchmark 
point has m ( - = 360 GeV and a representative DM mass 
= 50 GeV. We simulate t pair production with 
decay to a top quark pair plus dark matter candidates 
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at the LHC with 8 TeV energy. We demand that the 
top quark decays semi-leptonically, and that the anti-top 
quark decays hadronically, t — > 3jets. The final state 
contains a lepton plus jets and large missing transverse 
energy E*r- Two irreducible SM backgrounds, tt and 
ttZ production, are considered. Both the signal and 
background processes are generated at leading order in 
MadGraph/MadEvent Q with CTEQ6L1 parton distri- 
bution functions The renormalization and factoriza- 
tion scales are chosen as m^. Momentum smearing effects 
are included through a Gaussian-type energy resolution. 
We apply a set of basic acceptance cuts for the jets and 
single lepton in the final state: pt{C-) > 20 GeV, pr{j) > 
25 GeV, \r)t d \ < 2.5, A /■',.., , > 0.4, E T > 25 GeV. 
To suppress SM backgrounds, we impose a set of much 
harder cuts: prijist) > 50 GeV, pr{j2nd) > 40 GeV, 
E T > 100 GeV, H T > 500 GeV, where H T is the 
scalar sum of the transverse energies of all objects in 
the event. After the hard cuts, the cut efficiency for 
the signal is about 44% compared to the rate after the 
basic cuts. The ti background still dominates after the 
hard cuts, and the tiZ background is negligible. In order 
to further suppress the SM background, we use the fact 
that Et originates from the neutralino and neutrino in 
the signal events while from only the neutrino in the tt 
background. Hence, the neutrino longitudinal momen- 
tum pvL obtained from the VF-boson on-shell condition 
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is unphysical more often in the signal than in the back- 



ground 9]. Here A 



Hi' 



2p e T ■ Et- We then de- 



mand A 2 — ip 2 T E4 < 0. We also impose a cut on the 
transverse mass of the charged lepton and missing en- 
ergy, M T = \j2p e T Bt(1 - cos</>) > 100 GeV, where p T 
is the lepton transverse momentum and <f> is the angle 
in the transverse plane between px and E^- Only about 
0.00556% of the ti events remain after all the cuts. The 
cross sections for the signal and main backgrounds are 
shown in Table U after branching fractions are included. 
Using these cross sections, we find that the numbers of 
signal and background events are 130 and 22 at 8 TeV 
and 20 fb _1 integrated luminosity, for a signal signifi- 
cance of S/y/B = 28. 

In t pair production the decay chains of t — > t\ and 
t — > t\ have similar kinematics because the heavy i's 
are not highly boosted. In this work we investigate the 
energy of the anti top-quark as an estimator of the top 
quark energy, with the anti-top quark required to decay 
into three jets 



13j. We define a new energy fraction 



variable x' e , 

x' e = 2Et/E- t . (7) 
After convolution with the production cross section, a 



TABLE I: Cross sections (in fb) for the signal and back- 
grounds processes at different cut levels, including the decay 
branching fractions to the specific final states of interest. 





Basic 


thad recon. 


Hard 


Et sol. 


£cut 


signal 


22.26 


18.46 


8.87 


6.51 


11.6 % 


ti 


4347.08 


3596.75 


154.47 


0.91 


0.00556% 


tiZ 


1.25 


1.03 


0.34 


0.22 


5.9 % 



ratio 1Z' can be defined as 

1 /■*« da 
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cr(tot)7 
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dxf Xl ~ a(tot) ' (8j 



where da/dx' e is the differential cross section, and x' c is 
the cut threshold of the energy fraction x'g. 

We use a x 2 -template method based on the W boson 
and top quark masses to select the three jets from the 
hadronic decay of the anti-top quark. For each event we 
pick the combination which minimizes the following \ 2 '- 
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where Atovk an d Am t are the width of the lU-boson and 
the top quark, respectively. The efficiency of this method 
is 84%. After the antitop quark energy is reconstructed in 
the lab frame, 1Z' can be obtained with its cut threshold 
x' c dependence. 

Armed with both the Monte Carlo level momenta and 
the reconstructed momenta, we perform several compar- 
isons to evaluate how faithful the 1Z' distribution is to the 
true 1Z. At the Monte Carlo level, ti cp and fhad are known 
in the center-of-mass (cms) and lab frames. Our compar- 
isons show that 1Z defined with ii ep is not sensitive to the 
boost from the cms to laboratory frame, whereas 1Z de- 
fined by fhad shows a slight dependence. We compute the 
ratio 1Z defined from the energy of the t\ ep and thad- At 
the detector simulation level, only the four-momentum 
of fhad can be reconstructed, denoted i h ° a c d . Some of our 
results are compared in Fig. [5] (a) for choices sm^g = 1 
and cos0 c ff = 1 in Eq. O With sin# c ff = 1 (cos# c ff = 1) 
the top quark is mainly right-handed (left-handed), and 
we label the curves by tn (ti). There is some differ- 
ence between the 1Z distributions for t\ ep and fhad, but 
the essential features are preserved. We conclude that 
x' c is a good variable when x c cannot be obtained. We 
also investigate the cut dependence of f h ° a c d at the recon- 
struction level, whether basic or hard, and find that 1Z is 
not sensitive to the cuts; the curves for the loose cuts and 
the hard cuts overlap. Lastly, comparing 1Z at the Monte 
Carlo level and at the reconstruction level, we see a slight 
downward shift for both tj, and t^. This effect arises be- 
cause the pt cuts on the lepton reduce the number of 
events with x' t < x' c . 

The results in Fig. [5] (a) establish that x' c is a suitable 
variable and that 1Z' serves as a good substitute for 7Z. 
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FIG. 2: (a) The TZ distributions as a function of the cut 
threshold x' c for a 350 GeV t quark with pure right-handed or 
left-handed couplings at the LHC with 8 TeV energy. The lep- 
ton energy fraction is evaluated in the lab frame from either 
*£ad ( x 't) or the top quark energy t\ ep (xi). (b) The statis- 
tical uncertainty band of TZ is shown for both tL and tn at 
the lcr confidence level for an assumed 20 ftY 1 of integrated 
luminosity. 



To show that the difference between the expectations for 
fx and tR can be observed, we include statistical uncer- 
tainties on the TZ variable. Expressing 



n'{x' c ) = 



N(x' t < x' c ) 
Y(tot) ' 



(10) 



where iV(tot) is the total event number after cuts, and 
N(x' e < x' c ) (N(x' e > x' c )) represents events with x' e < x' c 
{x'p > x' c ), we derive the standard statistical uncertainty 



of 1Z can be used to identify the top quark polarization. 
In order to distinguish from an unpolarized top-quark, 
we define 10 bins in x' c within the range (0.3, 1.3) and cal- 
culate x 2 P er degree-of-freedom (d.o.f) for the difference 
n' L (x' c )-TZ' (x' c ): 
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The subscript "0" denotes an unpolarized top-quark. 
The result for an unpolarized top-quark (to) is 



K(x' c ) = {n' L {x' c ) + n'{x' c ))/2 



(13) 



After all the cuts, about 87 (101) signal events are needed 
to distinguish t L (t#) from t at 95% confidence level 
(C.L.) in the absence of background, with only about 23 
signal events to discriminate i# from t L . 

Thus far in this section, we assume the t-t-x coupling 
is completely left-handed or right-handed, but in general 
the coupling is a mixture of both. Once data are ob- 
tained, we could use the TZ' curves shown in Fig. 0(b) as 
templates in fits to these data to extract 6 c g and shed 
light on the nature of top squark mixing. 
Other implications: Our method can be applied to 
several NP models. We performed a detailed simula- 
tion of pair production of a T-odd top quark partner 
(T_) in the Littlest Higgs Model with T-parity (LHT), 
pp -> TJ?-X -> ttA H A H X, where A H is the T-odd 
photon partner. Our numerical results are very similar 
to those shown for tR in Fig. [21 Verification of mainly 
right-handed polarization would provide a powerful check 
of the model [ll|. Another example is the leptophobic 
Z' boson, which couples only to the SM quarks. The top 
quark polarization could be used to probe the handedness 
of the Z'-q-q coupling which is sensitive to how the SM 
quarks are gauged under the new gauge symmetry [12| . 
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The statistical uncertainties of the TZ distributions for tL 
and tR are shown in Fig.[2](b) at the lcr confidence level 
with 20 fb _1 integrated luminosity. The distinction is 
evident between L and R. The t L curve increases much 
faster than the tn curve. The pattern of the distributions 
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